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Abstract 

We study the effect of the late decaying saxino (the scalar superpartner of 
the axion) and find out that there is a possible dark matter solution from a 
class of supersymmetric extensions of the invisible axion model. In this class 
of models, the saxino which decays into two axions acts as the late decaying 
particle which reconciles the cold dark matter model with high values of the 
Hubble constant. Recent observations of the Hubble constant are converg- 
ing to Hq = 70 — 80kmsec~^ Mpc~^, which would be inconsistent with the 
standard mixed dark matter model. This class of models provides a plausible 
framework for the alternative cold dark matter plus late decaying particle 
model, with the interesting possibility that both cold dark matter and the 

extra radiation consist of axion. 
95.35.+d, 14.80.Ly 98.80.Cq 
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Cosmological and astrophysical data accumulated in past several years made it possible 
to test the theories of large scale structure formation in our universe. The COBE-DMR 
discovery of cosmic microwave background (CMB) anisotropy provided the means for the 
accurate normalization of theories of structure formation. The study of galaxy correlations 
supplied complementary data on smaller scales. With the scale invariant initial spectrum, 
they manifested that the pure hot dark matter model cannot explain the small scale structure 
of the universe and the pure cold dark matter(CDM) model shows best fit at the values 
floh = 0.25 [m {Qq =(the present energy density) /(the critical energy density) and h =(the 
present value of the Hubble constant) /lOO km sec~^ Mpc^^), which is much smaller than the 
presently favored values. The mixed dark matter(MDM) model attracted broad attention 
because it gives good fit for Q^, = 0.2 and ^cdm = 0.8 at h = 0.5 0]. 

For better test of models, however, the accurate values of Qq and h have been required. 
The observed value of Qq is ranging from 0.2 to ~ 1, still plagued by large systematic 
uncertainty. Several recent investigations of the Hubble constant tend to give higher values, 
in the range of h = 0.7 — 0.8 |^. If we stick to the theoretical prejudice, Qq = 1, we have 
at least two serious problems. First, we cannot fit the power spectrum curve even in the 
MDM model. Second, it would be inconsistent with the estimated lower bound on the age 
of the universe from the oldest globular clusters. An fio = 1 universe has an age of only 
to = Q.5h^^ Gyr, giving to < 9.3 Gyr for h > 0.7. On the other hand, the observed value of 
the age of the oldest globular cluster is around 15 Gyr |Q . Still there are a few possible way 
to relax this bound, but it seems not to be less than 11 Gyr ||5[. Assuming that the data 
from the globular cluster can be relaxed by some mechanism and considering error bars in 
Hubble constant observations, h = 0.6 would be marginally allowed. However, the standard 
MDM model is inconsistent with large scale structure data even for h = 0.6. 

Introduction of the small cosmological constant corresponding to = 0.7—0.8 alleviates 
the universe age crisis mildly and can revive good features of the CDM model with Q\ + 
^CDM = 1 (the CDM+A model). However, such a small value of cosmological constant is 
still a theoretically knotty subject. One way to keep the cosmological constant to be zero and 
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maintain good features of the CDM model is introducing the late decaying particle(LDP) 
(the CDM+LDP model) 0. In the CDM+LDP model, the LDP decays into very light and 
weakly interacting particles in the late stage of universe. This new radiation dominates the 
radiation energy of universe and would delay the beginning of the matter dominated epoch, 
which is necessary for the CDM model to be consistent with structure formation data even at 
the high values of h. But there are severe restrictions on the mass, lifetime and interactions 
of the LDP coming from structure formation data and preserving the successful predictions 
of nucleosynthesis and the CMB spectrum. We find that these can be met for the axion 
supermultiplet in supersymmetric extensions of the well-known axion model. It can provide 
cold dark matter with the late decaying particle consistently in a class of models. 

The axion is originally introduced to solve the strong CP problem. The axion model has 
an interesting cosmological effect. The axion produced from the initial vacuum misalignment 
has very small momentum and therefore is a good candidate of the cold dark matter. 

In the supersymmetric extension of the axion model, the axion supermultiplet ($) con- 
sists of the axion(a), its real scalar superpartner saxino(s), and the fermionic superpartner 
axino(a). The cosmological impacts of these two additional weakly interacting particles were 
studied before In this paper, we present unappreciated possibihty of saxino cosmology, 
in which the axino decays to two axions and acts as a LDP. 

After supersymmetry breaking, the axino and the saxino get their masses. In global 
supersymmetry, they remain massless at the tree-level. In the supergravity model with 
supersymmetry broken in the hidden sector, they are expected to gain masses of order 7713/2. 
However the axino and saxino masses are dependent on the specific forms of the axion sector 
superpotential as well as the Kahler function in the context of supergravity A model- 
dependent analysis is necessary to evaluate to the masses of axion supermultiplet. With the 
minimal Kahler function, the saxino mass is roughly of order of ^3/2 in most cases. In no 
scale supergravity, it remains massless at the tree level and gains radiative corrections of 
order 10 ~ lOOMeV. 

The effective Lagrangian describing the interactions of the axion supermultiplet is given 
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by Hi 



+ ^ {aF;;,F^'^^ + -a^.a.A^F^'^-' + sF^.^F'^^^ + ■■■)- -d,a ~a^^~a + ■ ■ ■ . (1) 

where Fa is the axion decay constant and etc is color couphng constant. The model-dependent 
parameter x is given by x = Y^illt^f I where and Vi are the [/(1)pq charges and the 
VEVs of the fields in the axion sector. It is of order 1, in general. At present particle 
phenomenology, astrophysical and cosmological observations restrict the range of the axion 
decay constant to be 

10^°GeV < F„ < lO^^GeV. (2) 

From these couplings we obtain the saxino decay widths to two different channels 

_ /o^ 
^^'^ 1287r3F2' ^""^ 

_ x^ml 

where is the saxino mass. Depending on the model dependent parameter x, the main 
decay mode can be changed and the effect of the saxino decay appears quite different. When 
s ^ 2g is the main decay mode, the lifetime is given by 

r. = 2.6 X 10 sec [-) [■^^^) [-^) ■ (5) 

The resulting gluons thermalize and dump some amount of entropy to thermal bath [§. 
Previous studies on the saxino decay have been focused on the saxino heavier than 10 GeV 
because it has been believed that the saxino should decay before nucleosynthesis. The reason 
is that the high energy radiation from the late decaying saxino would destroy the light nuclei 
which were made during the nucleosynthesis era. However when x ^ 10^^, as expected in 
many models, the main decay mode of the saxino is s — 2a with the lifetime 



r, = 1.3 X lO^sec x'^ , ° . (6) 

The crucial difference is that the resulting axions interact very weakly and cannot be ther- 
malized. They do not dump entropy to thermal bath of the universe and do no harm to 
the light nuclei made during nucleosynthesis. However, their effect appears gravitationally 
through their energy density. It affects the result of nucleosynthesis, which makes us consider 
the different range of the saxino mass. 

For further discussion, we need to know the relic abundance of the saxino. There are 
two distinct contributions: thermal and non-thermal. First we consider the thermal relic. 
Without inflation, the saxino would be hot thermal relic which decouples at the temperature, 
as estimated by Rajagopal et al. 0, 

Then its relic density is given by 



Ys = — = 0.278^^ - 1.2 X 10-^ (8) 

where Ug is the number density of the saxino and s = {27t'^ /45)g^sT^ is the entropy density. 

In the inflationary scenario, the reheating temperature is an important parameter in our 
discussion because the decoupling temperature of the saxino is rather high. At present, the 
relevant upper bound on the reheating temperature comes from the gravitino production 
and is quoted as [|rT| 



T,eh ^ lO^GeV (9) 

for reasonable values of the gravitino mass. If we use the MSSM value ac(lO^^GeV) ~ 1/20 
in the Eq. (|^), the reheating temperature seems lower than the decoupling temperature. 
However, in the heavy quark axion model, the heavy quark and saxino coupling is about 
(mass of heavy quark)/Fa, which is much stronger than any other interaction with saxino. 
In this case, the axion will not decouple until the heavy quark decouples from the heat 



bath. Therefore the decouphng temperature could be lower than reheating temperature. 
Though the decoupling temperature is greater than the reheating temperature, saxino could 
be produced from the thermal bath through the reactions like qq sg and gg ^ sg. Then 
the relic density is estimated to be 

y,.0.75xl0-(i;i|)(^)-'(^). (10) 

where 77 is phase space factor times the number of channels, and correction from the thermal 
effects which is roughly of order 1. 

As non-thermal relic, there can be a coherent oscillation of the saxino field. The coherent 
oscillation begins when the supersymmetry breaking saxino mass becomes comparable to 
the Hubble parameter. We estimate the rehc density from coherent oscillation to be 

where Si is the initial value of the saxino when the oscillation begins. The misalignment 
production is negligible for the saxino when Fa < 10^^ GeV. 

Prom the relic abundance of the saxino, we obtain strong bound on the saxino mass 
comes from nucleosynthesis. To preserve successful nucleosynthesis, the additional energy 
density at the time of nucleosynthesis should be smaller than, say, AA^,^ times the energy 
density of one neutrino species. For the saxinos and decay produced axions, this corresponds 
to mY < 1.3AN,{TD/g.siTn)) for To > T^s and mY < 1.3AN,{T^s/9*s{t^s)) fo^ To < T^s 



where = 0.55g^ToV MpT is the temperature at decay time and Tns — 0.8 MeV the 
temperature at the time of nucleosynthesis. Combining these, we obtain the axino mass 
bound 
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or 



< lOTMeVAAT, (^^) . (13) 



If we adopt AA^^, = 0.3 [|T^], the upper bound in the Eq. (|T3D would be 30MeV in the 



hot thermal relic case. Though the low reheating temperature weakens it, it is still difficult 
to get such a small mass in models where the saxino gets its mass at tree-level. However 
in models like no scale supergravity models, where the axino and the saxino have zero tree 
level masses and get masses by radiative corrections, the 10 ~ 100 MeV mass of the saxino 
naturally occurs in connection with the ~ keV order axino mass. 

Then the mass, lifetime, decay products and relic abundance of the saxino nicely fit with 
the requirements of the LDP. The condition the late decaying particle should satisfy is ^ 

3/2 
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For the hot saxino thermal relic, using the Eq. (Wj, it becomes 



(14) 



190 / \-i / Fg \' ^ 

x2 UoMeVy UOiiGeVy ' ^ ^ 

for the h = 0.6 case. Considering the bounds on Fa and (the Eqs. (|) and (|^)), this 
can be met in models with x ~ 1 and Fa ~ 10^° GeV. If the reheating temperature is lower 
than decoupling temperature of saxino, we should use the Eq. (plQl). Then we obtain 

7_0 ( ( Fa f r,eh \' _ . . 

x2 UoOMeVy VlO"GeVy \ 10^ GeV J ' ' ^ ^ 

In this case, depending on the reheating temperature T^eh, the condition can be satisfied in 
models with small x and Fa — 10^^ GeV. 

The CDM+LDP model is distinguished from the CDM+A model by the existence of the 
intermediate matter domination era before the final matter domination era. This temporary 
matter domination occurs because the saxino dominates the universe before it decays into 
the axions. Due to the temporary matter domination, there is one more length scale 

Aeqi ^ 8 X 10-2 (^) kpc. (17) 

Objects on this and smaller scales would collapse at high red shift and the CDM+LDP 
model has more structure on these scales. For hot saxino thermal relic, this scale is 



and for regenerated saxinos, 



The existence of this scale might be better in explaining the abundance of the high red shift 
quasars and globular clusters. 

We should check one more constraint on the LDP. Though the main decay mode of the 
saxino is s — 2a, the radiative branching ratio should be sufficiently suppressed not to 
spoil the nucleosynthesis result by photodestruction and photoproduction of light nuclei. 
For rris > 10 — 50MeV and > 10^ sec, it is roughly rrrisYs < lO^^MeV where r is the 
radiative branching ratio |13[. If the saxino is lighter than 100 MeV, the saxino cannot 



decay into hadrons, but can decay to two photons or two leptons. These decay channel is 
highly suppressed if the saxino is much heavier than 1 MeV. Decay to two photons is also 
suppressed by a^. Two lepton mode is suppressed by loop factor in the heavy quark axion 
model. 

Finally we discuss the possible CDM candidates in the above context. Actually any 
CDM is good. But because the supersymmetric axion model has two candidates, the axion 
and the axino, in itself, we discuss only about them. Though the mass of the axion is very 
small {rria — Aq^j-j/Fa), the large amount of axions can exist in our universe in the form 



of the coherent oscillation arising from the initial misalignment of the vacuum [jT4[. Its 
present energy density grows as Fa gets larger and becomes as large as the critical energy 
density around Fa = 10^^ GeV. The upper bound on Fa comes from the critical energy 
density. When the reheating temperature is larger than the decoupling temperature of the 
saxino. Fa ~ 10^'^ GeV is necessary for the saxino to be the LDP. But when the reheating 
temperature is lower. Fa — 10^^ GeV is desirable because the saxino can be the LDP even 
in models with small x. Therefore, we have an interesting possibility that CDM consists of 
the axions, and the saxino acts as the LDP, which decays to provide additional radiation 
which consists again of axions. 
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The decoupling temperature and the rehc number density of the axino are quite similar 
to those of the saxino. When Tj-eh > T^ec, the stable ~ 2keV mass axino can be so called 
warm dark matter When Treh is lower, the corresponding mass is 

and the axino becomes CDM. Interestingly enough, this axino mass range for which the 
axino is CDM and the saxino mass range for which the saxino is a LDP simultaneously 
occur in models where the axino and the saxino get their masses by radiative corrections 

m- 

In conclusion, we point out that the axion supermultiplet which solves the strong CP 
problem can provide cold dark matter with ^cdm = 1 and the late decaying particle simulta- 
neously within the presently allowed parameter range if the 10 MeV ~ 1 GeV saxino decays 
mainly into two axions. All the discussions are based on the flat, cosmological constant free 
universe. In this context, this model is one of a few possible models which satisfy the large 
scale structure and the large Hubble constant. It can survive the universe age problem if 
we allow the marginal value 11 Gyrs. If the cosmological and astrophysical observations are 
refined, the validity of this model will be proven in the near future. 
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